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ABSTRACT 

CMB experiments commonly use maps of Ha intensity as a spatial template for Galactic free-free 
emission, assuming a power law I v oc j/~ - 15 for the spectrum. Any departure from the assumed 
free-free spectrum could have a detrimental effect on determination of the primary CMB anisotropy. 
We show that the Ha-correlated emission spectrum in the diffuse WIM is not the expected free-free 
spectrum at WMAP frequencies. Instead, there is a broad bump in the spectrum at ~ 50 GHz which 
is consistent with emission from spinning dust grains. Spectra from both the full sky and smaller 
regions of interest are well fit by a superposition of a free-free and "warm ionized medium" Draine & 
Lazarian (1998) spinning dust model, shifted in frequency. The spinning dust emission is ~ 5 times 
weaker than the free-free component at 50 GHz, with the null hypothesis that the Ha-correlated 
spectrum is pure free-free, ruled out at > 8a in all regions and > lOOer for the full sky fit. 

Subject headings: diffuse radiation — dust, extinction — ISM: clouds — radiation mechanisms: non- 
thermal — radio continuum: ISM 



1. INTRODUCTION 

During the course of measuring the cosmic microwave 
background (CMB) anisotropy, the Wilkinson Microwave 
Anisotropy Probe (WMAP) has produced the most 
detailed and sensitive maps of interstellar medium 
(ISM) emission between 20 and 100 GHz taken to 
date. The quality of the data has profoundly influ- 
enced our understanding of the four main ISM con- 
tinuum emission me c hanisms in this f r equency range 
(IBennett et~aTl 120031: iFinkbeinerl l2003t iHinshaw et all 
120071 iDobler fc Finkbeinerl l2007t ). At the highest 
WMAP frequencies, the Galaxy is domi nated by thermal 
dust emission, traced by far IR maps (Finkbeine r et all 
1999). At lower frequencies, synchrotron radiation from 
supernova shock accelerated electrons spiraling in the 
Galactic magnetic field makes up the majority of the 
emission. In the inner Galaxy, the synchrotron emis- 
sion is brighter relative to low-frequency maps, and has 
a harder spectrum. 4 Next, free-free emission from un- 
bound electrons interacting with ions in ~ 10 4 K ionized 
gas has a harder spectrum than synchrotron, and is de- 
tectable in all bands. Finally, there is a bump in the spec- 
trum of dust-correlated e mission at ~ 20 GHz, sometime s 
called "Foreground X" (|de Oliveira-Costa et all |2002( ). 
Over the last decade, it has been recognized that this 
could be electric dipole emi ssion from rapidly rotatin g 
dust grains ( "spinning dust" ; IDraine fc Lazarianlll998bT) . 

1.1. Foreground X: spinning dust? 

Spinning dust is the electric dipole emission from the 
smallest dust grains which have a non-zero electric dipole 
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4 Finkbeiner (2004) has previously interpreted this as a separate 
componen t, the synchrotron "haze", but this interpreta tion is con- 
troversial {Finkbeiner 2004; Dobler & Finkbeiner 2007j). This haze 
component is spatially orthogonal to the other templates by con- 
struction, and though we include it in our fits below, its existence 
has negligible effect on the Ha-correlated spectrum. 



moment and are spun up by several different interactions 
with the interstellar medium (e.g., ion collisions, the 
interst ellar radiation field, etc.; see IDraine fc LazarianI 
Il998bf) . Spinning dust was considered as a potential 
sourc e of radio and microwave emission in o ther con- 
texts (lEricksonlll957tlFerrara fc Dettmarlll994[) . Because 
the model depends on the grain size distribution, ion- 
ization fraction, electric dipole moments, etc., the peak 
frequency can be anywhere from 10 — 50 GHz, but is 
usually in the 20 — 40 GHz range. This dependence of 
the spectrum on the ISM parameters has made spinning 
dust difficult to definitively identify, but may provide the 
qualitative behaviors necessary to make a fool-proof de- 
tection. 

In the COBE/DMR data, iKogut et alj (| 19961 ) found 
that the dust-correlated emission does not simply fall off 
monotonically with decreasing frequency as would be ex- 
pected for a thermal dust tail, but instead was greater at 
31 GHz than at 53 GHz. This be havior was confirmed by 



the S askatoon CMB experiment (|de Oliveira-Costa et al 



1997 
1997 



the 14.5 and 32 GHz OVRO data (Leitch et al 



the Cottingham fc B oughn 19.2 GHz survey 



(Ide Oliveira-Costa et al.lfl998l). t he Tenerife 10 & 15 GHz 
survey (|de Oliveira-Costa et alj I1999D. and the QMAP 
30 GHz data (|de Oliveira-Costa et all 120001 ) . A small 
survey of dust clouds with the Green Bank 4 2m tele- 
scope at 5, 8, & 10 GHz (jFinkbeiner et al.ll2002D revealed 
two candid ates for spinning dust clouds: L1622, later 
confirmed (jFinkbeinerl 120041) : an d LPH 201.663+1.643 , 
which turned out to be incorrect (Dickins on et al.| [2006) . 

Given this evidence, it was already clear prior to 
the WMAP data release in early 2003, that there 
was some non-standard emission mechanism, but it 
was unclear whether this m echanism was spinning dust 
(jDraine fc Lazarianlll998bl) or so me other process. 

WMAP ([Bennett et aTT 120031) found that the dust- 
correlated emission falls from 94 to 61 GHz and then 
rises steadily to 23 GHz. The fact that sy nchrotron typi- 
cally dominates at lower frequencies led IBennett et al.l 
(2003) to interpret the low frequency rise as "dust- 
correlated synchrotron" . The reasoning was that be- 
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cause dusty regions correlate with regions of star for- 
mation activity, electrons in the vicinity could be accel- 
erate d by supernova sh ocks leading to synchrotron emis- 
sion. iFinkbeinerl (120041) pointed out that W MAP is also 
consistent with lDraine fc Lazarianl (|1998bh if their model 
could be shifted down in frequency a bit. 

Because the WMAP frequency coverage (23 — 94 
GHz) coincides with the range where the spinning 
dust and synchrotron spectra are similar, data at 
lower frequencies, where the spectra diverge, have been 
used to test the dust-correlated synchrotron hypothe- 
sis. The WMAP "synchrotron" template has a much 
weaker cross-correlation at 10 — 15 GHz than expecte d 
for actual synchrotron (|de Oliveira-Costa et al.l |2004[ ). 
Also, 8 & 14 GH z data from the Gree n Bank Galac- 
tic Plane Survey (|Langston et alj l2000f ) show a rising 
spectrum consistent with spinning dust but not syn- 
chrotron (|Finkbeiner. Langston. fe M intcr 2004j) and im- 
plying that the majority of dust-correlated emission at 
23 GHz is Foreground X. Comparing the Cottingham & 
Boughn survey to WMAP also shows a rise from 19 to 
23 GH z dBoughn fc Poberll2007l). The polari zation prop- 
erties (page et al.ll2007t iHinshaw et al.ll2007l ) in the 3-yr 
WMAP data also favor the spinning dust hypothesis. 

In spite of all this evidence, there are weaknesses 
in each of these arguments. Some use only low- 
latitude data and may not apply to the diffuse ISM 
(|Finkbeiner. Langston. fc Mintei I2004D . Others apply to 
only one cloud (Finkbci ner et alJ 20021) , have poor signal 
to noise (|de Oliveira-Costa et al.ll2004h. or a short lever 
arm in frequency l|Boughn fc Poberl |2007| ). However, 
taken together, they are very strong evidence against the 
dust-correlated synchrotron hypothesis, but they are not 
specifically evidence for spinning dust. Other mecha- 
nisms such as magnetic dipole dust emission, caused by 
the the rmal fluctuation of the m agnetization in the dust 
grains (jDraine fc Lazar ian 1999), could conceivably ex- 
plain all these results. This ambiguity motivated our 
search for further clues in the WMAP 3-yr data. 

1.2. WMAP Multilinear Regression 

In lDobler fc Finkbeinerl (|2007l . hereafter, DF07) we de- 
scribe our method for determining the spectral shape 
of each of the four primary WMAP foreground compo- 
nents: free-free, dust, soft sync hrotron, and har d syn- 
chrotron (termed the "haze" by iFinkbein cr 2004) emis- 
sion. We model the WMAP data as a linear combination 
of the CMB and the 4 foregrounds, each described by a 
specific spatial template. For free-free emission we use 
the Ha map described in Q2.1\ for dust (therma l and 
spinning) emission we use the|Schlcgcl et al. (1998 1) dus t 
map evaluated at 94 GHz bv [Finkbeiner et al.l (11 999). 
for soft synchrotron we use the lHaslam et al.l (|1982h 408 
MHz map, and for the haze we use a simple 1/r tem- 
plate where r is the distance from the Galactic center. 
The individual WMAP bands were completely decou- 
pled in our fit, and we assumed nothing about the spec- 
trum of each foreground, only that its morphology was 
traced by the template in all bands. We performed our 
fits over both the (nearly) full sky and also on smaller 
regions of interest. Point sources and regions where 
dust extinction makes the Ha map unreliable (where 
A (Ha) = 2.65 E(B — V)>1 mag) were masked out from 
the fits. 



We showed that the spectrum of Ha-correlated emis- 
sion is not the usual free-free spectrum. Rather, there 
is a bump around v ~ 50 GHz. In ij2.3l we explain why 
the Ha map should actually trace spinning dust excited 
by collisions with ions. In Sj3] we present Ha-correlated 
emission spectra for both the (nearly) full sky as well 
as smaller regions of interest and show that, indeed, the 
spectra are well characterized as a superposition of free- 
free and spinning dust type spectra. We summarize this 
work and discuss implications for future observations in 

m 

2. EMISSION MEASURE AS A TRACER OF SPINNING 
DUST 

In this section we argue that in fully ionized environ- 
ments, the emissivity density (j u ) of spinning dust emis- 
sion is proportional to dust density times some excitation 
function related to the ion density. The emission comes 
from the very smallest grains, which are spun up by a 
collision with an ion and spin down as they radiate. In 
the limit where the spin down time is small compared to 
the time between collisions, the rotation is episodic and 
the emissivity also scales linearly with the ion number 
density. In ij'2.21 we assume the episodic limit in order to 
build some intuition about how WIM-correlated spinning 
dust might behave. 

2.1. Free-free emission measure 
The map assembled bv | Finkbeinerl (l2003f) from the 



VTS S (jDennison et all 11998ft. SH ASSA (fGaustad et all 
120011 ). and WHAM (jHaffner et al.ll2003ft surveys, traces 
the Ha recombination line emission from hot interstel- 
lar gas. At WMAP frequencies, this Ha map 5 is used 
to estimate the thermal bremsstrahlung, or free-free 
emission, resulting from the interac tion of the electrons 
and proton s from ionized hydrogen dBennett et aTll2003l : 
IFinkbeinerl 12004 IHinshaw et all 120071 ). The Ha map is 
corrected for dust absorption using the prescription given 
by Finkbeiner (2003; but see Dickinson et al. 2003 for an 
alternative) . In this frequency range the free-free specific 
intensity l v (in Jy/sr) is often approximated as a power 
law, 

(1) 

with a « —0.15. Although this is a good approxi mation, 
in our analysis we use the expression given by iSpitzerl 
(|1978ft for the free-free spectrum. 

Recombination line emission is a two particle process 
(involving the capture of free electrons by a proton nu- 
cleus) and so the emissivity is proportional to n e rii, which 
for fully ionized H is simply the square of the number 
density of gas particles, Jjjq, oc n 2 . Therefore, the Ha 
intensity corresponds to an emission measure (EM) de- 
fined as E m = J n 2 dl where the integral is taken along 
the line of sight. The Ha intensity is usually measured in 
Rayleighs (1 R = 10 6 /47r photons s _1 cm~ 2 sr _1 ) with 
an EM of 1 cm- 6 pc ps 0.6 R for gas at - 8000 K. 

2.2. Spinning dust emission measure 

iDraine fc Lazarianl (|1998bl hereafter, DL98) present 
a model for spinning dust emission in which tiny dust 
grains with non-negligible electric dipole moment are 

5 available at http://www.skymaps.info 



3 



spun up by ion collisions; polycyclic aromatic hydrocar- 
bons (PAHs) generally have a non-zero dipole, because 
perfect symmetry is rare in large molecules. DL98 esti- 
mated the astronomical value of this parameter (actually 
the electric dipole moment per root mass) from mea- 
sured electric dipole moments for dozens of known PAH 
molecules. In the limit where ion interactions provide the 
dominant spin-up mechanism (see DL98 for details) and 
assuming full ionization, constant gas/dust ratio, and the 
same grain size distribution everywhere, the resulting to- 
tal spinning dust emission vs I = J l v dv oc J Uin d dl oc 
J n 2 . With the additional assumption that the emission 
is episodic, the spectrum of the emission does not vary 
with m, but the amplitude does. That is, the spinning 
dust emission should look like an EM, and be traced by 
Ha. This approximation is relevant where the collisional 
excitation is primarily due to ions (at least for the small- 
est grains), e.g. for the "warm ionized medium" (WIM) 
parameters used by DL98. 

The WIM-correlated emission claimed here augments 
the "usual" spinning dust emission in the "cold neu- 
tral medium" (CNM) which has a temperature Tq nm ~ 
100 K and is traced by the iSchlegel et "all (119981) dust 
map evaluated at 94 GHz by iFinkbeiner et al.l (|1999f ) 
(FSD99). 

2.3. Estimating the specific intensity 

For a spherical dust grain of size a, rotating with an- 
gular frequency u>, and with a dipole moment po, the 
rotational energy is, 



E = -Ioj 2 
2 



1 2 2 

—ma lu , 



where m is the grain mass. The energy loss rate is the 
total power emitted by the grain, 



TP- 1 ^°„2 4 

12ir c 



(3) 



and the spin down time is, 

t s = E/E oc a 2 P Q 2 uj- 2 . (4) 

The time between collisions is just, 

r c ~ A/u a a~ 2 n~j^ (5) 

where v is the thermal velocity of the ion and A = 
(7ra 2 nj^) _1 is the mean free path. 

In the limit where the spin down time is much smaller 
than the time between collisions (r = t s /t c <C 1), the 
excitation is episodic. For the WIM parameters in DL98, 
this ratio is in fact ~ 1; however, we point out that 

r = t s /t c oc a 4 n H pQ 2 uj~ 2 (6) 

so that very small changes to the grain size distribution 
can substantially affect the spin down time compared to 
the time between collisions. In fact, if each ion collision 
imparts equal energy, then from Equation [2l lo oc a~ 2 . 
Since po oc TV 1 / 2 oc a 3 / 2 (see DL98) then in this case, 1/r 
diverges as a~ 5 . The bottom line is that the episodic 
approximation is sensitively dependent on the ion den- 
sity, the dipole moment, and especially the grain size 
distribution. For the purposes of estimating a specific 
intensity, we will now assume the episodic approximation 
with the caveat that it may break down in dense env iron- 
ments, (e.g., HII regions. In fact, iScaife et al.ll2007L find 



no evidence for spinning dust at 14.2 and 17.9 GHZ in 16 
Galactic HII regions) or for very small dipole moments. 
We emphasize that this does not affect the fact the the 
total intensity I should be proportional to J n d njid£ so 
that spinning dust emission should be traced by an EM 
map. 

Draine & Lazarian evaluate j v /nH for parameters rep- 
resentative of various interstellar environments. In the 
limit that j^/n^ is a constant with n#, then it equals 
I V /N(H). However, in our case it is not. As a refer- 
ence model, we take the DL98 WIM model evaluated 
for uh = 0.1 cm -3 which peaks near 25 GHz with a 
peak emissivity 6 of jp^cf/nn = 5 x 10 -21 kJysr~ 1 cm 2 or 
jp.rd = 5x 10~ 22 kJysr~ cm -1 . Multiplying the depen- 
dence on dust and gas density we have 7 



Jp — Jp^ci 



nH,ref 



n d 

n d,ref 



(7) 



where nH.rtf and n^rcf are the density of H and small 
dust grains, respectively, used in the reference model. 
For full ionization of H (jih — n^) and constant gas/dust 
ratio {n H /n H ^ c{ = n d /n d ^ ci ) we have 



• - ■ ( Uh V 
3p ~ Jp ' Iei \0.1cm- 3 J 

= 5 x 10~ 20 7i^ kJysr^cm 5 
= 0.15n^ kjysr _1 cm 6 pc _1 
and the peak specific intensity is then 



(2) I P = J jpdl = 0.15 kJysr- 1 ^~ 



E„ 



cm °pc 



or 



V 7 Ha = °- 09 kJysr^IT 



(8) 

(9) 
(10) 

(11) 
(12) 

where I^{ a is measured in R. For comparison, the ratio 
of free-free to Ha for a temperature of T e ~ 8000K at 41 
GHz (WMAP Q-band) is hx/I^ a - 0.15 kJysr" 1 R -1 . 

Although this is a rough estimate, it suggests that the 
intensity of the free free emission at 41 GHz and the peak 
of the spinning dust emission should be of the same order 
of magnitude. Throughout this discussion we make the 
assumption that the DL98 WIM spectrum can be shifted 
by a factor of 2 in frequency. A more thorough analysis 
would require that the DL98 model be evaluated for some 
set of parameters (ion density, electric dipole moment per 
root mass, etc.) that give the actual spectrum observed. 
For our purposes this naive shift is adequate. 

3. SPECTRUM OF Ha-CORRELATED EMISSION 

As noted above, in DF07 we showed that the spec- 
trum of Ha-correlated emission in the WMAP data has 
a bump around v ~ 50 GHz that we now argue is consis- 
tent with a spinning dust contamination of the familiar 
Equation [T] free-free spectrum. We also point out in 
DF07 that any measurement of a foreground spectrum 
is necessarily contaminated by a component which has 

6 The units of the ordinate, jv/nH, of Figs. 9-11 in DL98 is 
incorrectly labeled (Jy sr _1 ) rather than (Jy cm 2 sr — 1 ), an error 
propagated by Finkbeiner et al. (2002). 

7 Again, the dependence on njj in Equation [7] is only strictly 
correct in the episodic limit, though in reality the deviation is likely 
small. 
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Fig. 1. — Top panel: The best F+ D + C (free-free plus spinning 
dust plus CMB) model fit to the Hct-correlated emission over the 
full sky. The peak frequency is fixed to !^ pca k = 50 GHz. This two 
component plus CMB contamination model works surprisingly well 
with a spinning dust component that is 6.23 times weaker than the 
free-free (see Table IT) . Middle panel: A\ 2 contours for the best 
fit free-free coefficient in the spinning dust/CMB plane. The null 
hypothesis that the Ha-correlated emission can be well fit by only 
a F + C model (i.e., Drj = 0) is ruled at very high confidence 
(335<t). Bottom panel: x 2 contours in the spinning dust / f pea k 
plane. Though the peak frequency is not very well constrained, 
we can rule out f pea k < 40 GHz and f pea k > 60 GHz at > 4<r 
confidence. 



the frequency dependence of the CMB due to imperfect 
foreground removal when constructing a CMB template. 8 
However, as we will show shortly, this contamination is 
incapable of producing the observed bump. 

Figure [T] shows our full sky fit for Ha-correlated emis- 
sion. Our hypothesis is that this spectrum can be well fit 
by a linear combination of free-free, spinning dust, and 
CMB spectra. Thus, our model spectrum is 



rmod = F 



D + C, 



where 



F 



»( 



23 GHz 



-0.15 



(13) 



(14) 



D = D a x (DL98 WIM model with 



and 



C = C, 



50 GHz), 

(15) 

.23GHZ7 " (16) 

We minimize x 2 = J2i(^° d ~ ^i) 2 /°f over the parame- 
ters F , D , and Co which are scaled to be of order unity. 
Here <7j are the formal errors in band i of the spectrum 
from our DF07 template fit, and we have one degree of 
freedom. 

As shown in Figure [TJ the spectral model in Equa- 
tion [13] fits the data astonishingly well. The free-free 
and spinning dust amplitudes are comparable while the 
contamination from the CMB is minimal. The very low 
value of x 2 = 1-84 reflects the quality of the fit. Fig- 
ure [T] also shows A\ 2 = 1,2,3,4, and 5 contours in the 
Do-Co plane minimizing \ 2 over Fq at each point, and 
holding v p fixed at 50 GHz. The null hypothesis that 
only a free-free plus CMB spectrum (i.e., D = 0) fit 
the data is ruled out to a very high significance. Thus, 
either free-free emission does not follow Equation [TJ or 
the Ha map is tracing an additional emission component. 
This other emission component seems to be well fit by a 
shifted WIM spinning dust spectrum. 

Lastly, the x 2 contours in the -Do-^peak plane in the 
bottom panel of Figure [T] indicate that, although the 
peak frequency of the spinning dust is not very well con- 
strained, values of ~ 45 — 55 GHz are consistent with 
the data. Peak frequencies less than 40 GHz and greater 
than 60 GHz are ruled out at 4cr. 

The next step is to consider smaller individual regions 
of the sky to show that the signal persists. Figures[2]and 
[3] show fits for four different regions of the Gum Nebula, 
the north central Galactic sky (including the Ophiuchus 
complex) , and a region near the Orion Nebula (see Table 
[1] and Figure 2] for details). Though the fits are signifi- 
cantly noisier as a result of the smaller number of pixels 
used, the bump in the Ha-correlated emission remains 
in all regions, with the free-free amplitude roughly 4-7 
times larger than the spinning dust. Furthermore, the 
case of no spinning dust D = is ruled out in all re- 
gions at > 8a. Although the x 2 contours in the -Dp-^poak 
plane offer a tantalizing hint that the peak frequency of 
the spinning dust is lower in the Orion Nebula region 
than the others, the present data are simply too noisy to 
make a statistically significant statement. 

Lastly, we address the possibility that the bump in 
the Ha-correlated emission represents a cross-correlation 

8 Throughout this paper, we will work with the results found 
using our CMB5 estimator for the CMB - see DF07 for details. 
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Fig. 2. — The same as Figure [T] except for the four regions of the Gum Nebula (see Table [TJ. Although the peak frequency is only 
loosely constrained, the Do = null hypothesis is ruled out at > 12<r in all regions. 



TABLE 1 





Full sky 


1 


Gum I 
2 


Nebula 
3 


4 


Ophiuchus 


Orion 


lrange 


0:360 


235:256 


256:280 


235:256 


256:280 


334:13 


195:230 


brange 


-90:90 


-30:0 


-30:0 


0:30 


0:30 


0:90 


-35:0 


x 1 


1.84 


0.17 


1.10 


1.08 


0.26 


0.17 


1.52 




335.02 


43.84 


87.32 


17.59 


12.09 


11.89 


8.30 




6.23 


4.70 


6.07 


5.38 


6.80 


4.75 


12.22 



Note. — Region definitions (in Galactic coordinates I and ft) and results of the fits presented 
in Figures IU3I The very low values of \ 2 reflect both the quality of the fits and the small 
number of constraints (number of degrees of freedom = 2). The Do = null hypothesis is ruled 
out at > 8cr in all regions, while the relative intensity of free— free to spinning dust emission at 
^peak = 50 GHz is in the range 4.7-12.2. 



TABLE 2 





Ha 


Dust 


Haslam 


Haze 


Ha 


1.00 


0.45 


0.18 


0.10 


Dust 


0.45 


1.00 


0.56 


0.26 


Haslam 


0.18 


0.56 


1.00 


0.58 


Haze 


0.10 


0.26 


0.58 


1.00 


Note. 




The 


correlation ma- 



bias between the different templates used in DF07. That 
is to say, perhaps the Ha spectrum is absorbing power 
from other templates in the fit. Table [2] shows the corre- 
lation matrix, defined as 



(TiTj) 



{T?)(T?) 



(17) 



trix for the temp l ates presented in 
IDobler fe Finkbeinerl | |2007h . These are 
the templates that were used in the 
multi-linear regression fit which revealed 
the bump in the Ha-co rrelated emission 
presented in Figures |H3I 



where Fij is the correlation between the ith and jth tem- 
plate. The correlation between the templates is rather 
small. In particular, the dust map is the template which 
is most highly correlated with the Ha map, and even 
then, the variance in the Ha-correlated emission can be 
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Fig. 3. — The same as Figure[T]except for a region which encloses 
the Ophiuchus Complex (left column) and one which is near the 
Orion Nebula (right column). See Table ITl for details. 



at most r?j , , 
Ha, dust 



20% due to correlation with dust. 

In fact, Figure [5] shows that the effects are much 
smaller. Here, the solid lines show the resultant spec- 
tra for all foreground components both with and with- 
out fitting the Ha spectrum. In the latter case the Ha 
emission was assumed to follow the free-free spectrum in 
Equation [T] The features of these spectra are discussed 
at length in DF07, but the important point is that the 
other foregrounds are not significantly affected when the 
Ha spectrum is fixed. This is strong evidence that the 
emission mechanism causing the bump in the Ha spec- 
trum is highly correlated with the Ha template. 

4. DISCUSSION 

In iDobler fe Finkbeineri (|2007f ) we identified a 
deviation from the classical free-free (thermal 
bremsstrahlung) spectrum in the Ha-correlated emission 
of the 3-year WMAP data. In this paper we have argued 
that the spectral data are consistent with a superpo- 
sition of free-free and spinning du st emission from a 
"warm ionized medium" (WIM; see iDraine fc Lazariar] 
1998b). Additionally, there is contamination by a 
component with the spectrum of the CMB due to 
chance morphologi cal correlation between th e CMB and 
the Ha map f see IDobler fc Finkb cincr 2007], for details). 



Since the intensity of both free-free and WIM spin- 
ning dust emission scale proportionately with the num- 
ber density squared (free-free emission is generated by 
the collision of free electrons with ions while the WIM 
spinning dust emission comes from collision of ions with 
tiny dust grains), both should be traced by an emission 
measure map like the Ha map. 

We have assumed throughout that the dust/gas ratio 
is constant and the grain size distribution is the same 
everywhere, so that / oc J nin^dl oc J nf. 

We have studied both the (nearly) full sky and smaller 
independent regions of interest including the Gum Neb- 
ula, the Oph complex, and the Orion Nebula. For all 
regions, the Ha-correlated emission is well fit by a linear 
c ombination of an l v o c i/~ 01 5 power law for free-free and 
a IDraine fc Lazarianl (|T998b) WIM spinning dust model 
shifted in frequency by a factor of 2 and with only a min- 
imal adjustment in amplitude. Specifically, the spinning 
dust component is ~ 4-7 times weaker than the free-free 
component. 

The null hypothesis that a free-free only spectrum is 
consistent with the data is ruled out at > 8<r confi- 
dence for all regions considered here, despite the sub- 
stantial noise in the data. Furthermore, we have shown 
that the bump in the Ha-correlated emission spectrum 
is not due to chance cross correlation between the Ha 
m ap and the additional foregr ound templates presented 
in IDobler fc Finkbeineri (|2007f ). indicating that it is in- 
d eed correlated primarily with Ha emission. 

Finkbci ner et al.l (|2002l ) listed a number of criteria by 
which the existence of spinning dust could be established. 
Among them were that the peak frequency must shift as 
a function of the environment in a way at least qualita- 
tively consistent with DL98. In this work we find in the 
WMAP data a WIM-correlated component with a peak 
frequency of 50 GHz, about a factor of two higher than 
the peak frequency of the spinning dust component that 
is spatially correlated with thermal dust emission (i.e. 
WNM or CNM spinning dust). That this frequency is 
higher is in accordance with the DL98 models, strongly 
supporting the spinning dust hypothesis. However, be- 
cause the shift is larger than expected, a careful assess- 
ment of the situation requires that the reason for the 
size of the shift be determined. If, for example, the WIM 
dust has smaller grains, then a correlated measurement 
of Ry = Av/Eb-v using e.g. background stars would 
complete the case for spinning dust. 

In summary, there is clearly an additional foreground 
component correlated with Ha emission and its spectral 
shape is consistent with spinning dust, though at a some- 
what higher frequency than expected. However, the in- 
creased sensitivity, frequency coverage, and finer spatial 
resolution of future missions like Planck may provide the 
necessary evidence. 

We acknowledge helpful discussions with Bruce Draine, 
Gary Hinshaw, Joanna Dunkley, et al. Some of the 
results in th i s pap er were derived using HEALPix 9 
(jGorski et alJ 120051 ). This research made use of the 
NASA Astrophysics Data System (ADS) and the IDL As- 
tronomy User's Library at Goddard 10 . DPF and GD are 
supported in part by NASA LTSA grant NAG5-12972. 
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Fig. 4. — WMAP K band data (top) and the Ho map (bottom) with our different fit regions shown as dotted boxes. In addition to the 
full sky fit, the regions are the Gum Nebula (1,2,3,4), the Ophiuchus complex (5), and the Orion Nebula (6). See Table IT1 for details. 



9 The HEALPix web page is http://healpix.jpl.nasa.gov 

10 Available at http://idlastro.gsfc.nasa.gov 
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Fig. 5. — The inten sity of the emission correlat ed with the four 
templates presented in Doblcr & Finkbcincr (2007). The solid lines 
indicate fits for which the Ha-correlated emission was constrained 
to follow a Iff oc v~ - 15 (see Equation [TJ frequency dependence, 
while the dotted lines indicate that the Hct-correlated spectrum was 
fit explicitly. Both the insensitivity of the inferred dust, Haslam, 
and haze -correlated spectra to fixing the Ha spectrum as well as 
the relatively low correlation among the templates (Tabled indi- 
cates that the spinning dust bump is indeed correlated primarily 
with Ha emission. 



